Semiclassical theory of spin-polarized shot noise in mesoscopic diffusive conductors 
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We study fluctuations of spin-polarized currents in a three-terminal spin- valve system consisting of a diffusive 
normal metal wire connected by tunnel junctions to three ferromagnetic terminals. Based on a spin-dependent 
Boltzmann-Langevin equation, we develop a semiclassical theory of charge and spin currents and the correla- 
tions of the currents fluctuations. In the three terminal system, we show that current fluctuations are strongly 
affected by the spin-flip scattering in the normal metal and the spin polarizations of the terminals, which may 
point in different directions. We analyze the dependence of the shot noise and the cross-correlations on the 
spin-flip scattering rate in the full range of the spin polarizations and for different magnetic configurations. Our 
result demonstrate that noise measurements in multi-terminal devices allow to determine the spin-flip scattering 
rate by changing the polarizations of ferromagnetic terminals. 



I. INTRODUCTION 

Spin-dependent electronic transport and use of the spin de- 
gree of freedom of electrons in hybrid magnetic structures 
have recently been subject of highly interesting field termed 
spintronics. This developing field has emerged from many 
exciting phenomena, the giant magneto-resistance being the 
most well known example, which has made it attractive for 
applications as well as fundamental studies 1 - 2 - 3,4,5,6 - 7 - 8 . Typ- 
ical studied systems are based on the use of ferromagnetic 
metals (F) and/or external magnetic fields to inject, manip- 
ulate and detect spin-polarized electrons inside a mesoscopic 
normal metal (N). 

An important characteristic of mesoscopic systems is the 
appearance of shot noise, the fluctuations of current through 
the system due to the randomness of the electronic scattering 
processes and the quantum statistics. Shot noise and nonlo- 
cal correlations of the current fluctuations contain additional 
information on the conduction process which is not gain- 
able through a mean current measurement. In ferromagnetic- 
normal metal structures, in which the spin of electrons plays 
an essential rule, current fluctuations are due to the ran- 
domness of both charge and spin transport processes. Thus, 
shot noise measurements are expected to provide information 
about spin-dependent scattering processes and spin accumu- 
lation in the system. Together with the importance of noise in 
spintronic devices in view of applications this motivates our 
study of spin-polarized shot noise. 

In the past years, shot noise has been extensively studied 
in a wide variety of hybrid structures involving normal met- 
als, semiconductors and superconductors 9,10 . However, there 
are few studies devoted to shot noise in ferromagnet-normal 
metal systems. Results of earlier studies of current fluctua- 
tions in FNF double barrier systems in the Coulomb blockade 
regime 11 and FIF (I being an insulator) systems^ can be un- 
derstood in terms of well known results for the corresponding 
normal-metal systems for two spin directions 9 . Tserkovnyak 
and Brataas studied shot noise in double barrier FNF sys- 
tems with noncollinear magnetizations in F- terminals 13 . They 
found that the shot noise has a non-monotonic behavior with 
respect to the relative angle of the magnetizations for different 



scattering regimes and different types of FN junctions. 

The effect of spin flip scattering on spin-polariz ed current 
fluctuations has been considered in Refs. Il4ll5ll6ll7ll1 
Mishchenkc 14 found that in a perfectly polarized two-terminal 
double barrier system spin-flip scattering leads to a strong de- 
pendence of shot noise on the relative orientation of the mag- 
netizations in F-terminals. In Ref. 18 we have proposed a 
four- terminal spin- valve system of tunnel junctions to study 
spin-dependent shot noise and cross- correlations simultane- 
ously. It has been found that the cross-correlations between 
currents in terminals with opposite spin polarization can be 
used to measure directly the spin-flip scattering rate. 

Recently, there have been also studies of the current fluctu- 
ations of spin-polarized entangled electrons in quantum dots 
and wires 19 and in quantum dots attached to the ferromag- 
netic leads in the Coulomb blockade 20 , Kondc 21 and sequen- 
tial tunnelingi 2 * 2 ^!* 2 ^ 6 . regimes. 

In this paper, we study current fluctuations in a three ter- 
minal diffusive FNF system in the full range of spin polar- 
izations and the spin-flip scattering intensity. Based on the 
Boltzmann-Langevin 27,28,29,30 kinetic approach, we develop a 
semiclassical theory for spin-polarized transport in the pres- 
ence of the spin-flip scattering. We obtain the basic equations 
of charge and spin transport, which allow the calculation of 
mean currents and the correlations of current fluctuations in 
multi-terminal diffusive systems. Applying the developed for- 
malism to a three-terminal geometry, we find that current cor- 
relations are affected strongly by spin-flip scattering and spin 
polarizations. We focus on the shot noise of the total current 
through the system and the cross-correlations measured be- 
tween currents of two terminals. We demonstrate how these 
correlations deviate from the noise characteristics of the un- 
polarized system, depending on the spin-flip scattering rate, 
the polarizations of the terminals and their magnetic configu- 
rations (relative directions). Our results provide a full analysis 
for spin-dependent shot noise and cross-correlations in terms 
of the relevant parameters. 

The structure of the paper is as follows. In section II we ex- 
tend the Boltzmann-Langevin equations to the diffusive sys- 
tems, in which spin-flip scattering takes place and which are 
connected to ferromagnetic terminals. We find the basic equa- 
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tions of the charge and the spin currents and the correlations 
of their fluctuations. In section III we apply this formalism 
to the three terminal system. We obtain the Fano factor and 
the cross-correlations between currents through two different 
terminals. Section VI is devoted to the analysis of the cal- 
culated quantities. We present analytical expressions for the 
Fano factor and the cross-correlations in different important 
limits. Finally, we end with some conclusion in section V. 



II. BOLTZMANN-LANGEVIN EQUATIONS WITH 
SPIN-FLIP SCATTERING 

In this section we extend the semiclassical Boltzmann- 
Langevin kinetic approach 9 to cover spin-polarized trans- 
port. In the presence of spin-flip scattering the Boltzmann- 
Langevin equation is written as 



dt 



fa = I lmP [fa] + %[f a , f—a] + C P + Cl 



(1) 



where the fluctuating distribution function of spin a(= ±1) 
electrons f a (p, r, t) = f a (p, r) + 6f a (p, r, t) depends on the 
momentum p, the position r, and the time t. Eq. Q contains 
both normal impurity and spin-flip collision integrals which 
are given by the relations 

r mp {f a ] = nj j0-[j aa (p> >p ) - j QQ (p, P ')], (2) 

#[/«,/-*] =SlJ j0^[J- aa ( P ',p) - J„_ a (p,p')]- 

(3) 

Here J aa > (p, p', r, t) = W aa > (p,p' 7 r)f a (p,r,t)[l ~ 
/a'(p') r, t)], where W aa > (p, p', r) is the elastic scattering 
rate from the state p, a into p', a'; ft is the volume of the sys- 
tem. The corresponding Langevin sources of fluctuations due 
to the random character of the electron scattering are given by 

C P = ^ / ^0^[SJ aa (p',p)-SJ aa (p,p')], (4) 

£ = n J j0T3[SJ-a a (p',p)-SJ a - a (p,p%(5) 

where the random variable SJ aa /(p',p ) T,t) is 
the fluctuation of the current J aa i (p, p', r, t) = 

Jaa' (P, P', r, t)_+ 5 J aa ' (p, p', T, t)) with J aa > (p, p', r, t) = 

W aa >(p,p',r)f a (p,r,t)[l - f a '(p',r,t)] being the mean 
current. 

We will assume that all scattering events are independent 
elementary processes and thus the correlator of the current 
fluctuations 5 J aa ' (p, p', r, t) is that of a Poissonian process: 

< 5Ja 1 a 2 {Pl,P2,r,t)SJa 3 a 4 (P3,P4,r',t') > = 
(27T/i) 6 



{> -S ai a 3 Sa2aJ{Pl ~ P3)<5(p2 - P4) 

xS(r - r')S(t - t')J ai a 2 (pi, P2, r, t). 



(6) 



Due to the non-vanishing spin-flip collision integral in Q, 
the equations for the distributions of electrons with opposite 



spin directions are coupled. The coupled equations can be 
transformed into two decoupled equations for the charge / c = 
Tliafa/^ and spin / s = J2 a a fa/^ distribution functions, 
which read 

|/c ( s) = / imp [/c (s) ] + 4)[/c( S) ] + e:; p } + 4). (7) 

Here we have introduced different collision integrals: 



^ f [/c] = sij -0^ [Jf(p',p) Jf(p,P% 

Ififs] 



-n 



^3^ Sf (p,p')[/ S (p') + /s(p)], 



(8) 
(9) 

(10) 



where 



j5«(p,p',r,f) = irfW(p,p',r) 



xfc( s )(p,r,t)[l - /c(,)(p',r,i)], 



(11) 



and we assumed W aa '(p, p') = W a ' a (p' ,p), W-\ = 

W-+ = W st ; W++ = W— = W im P. The corresponding 
Langevin sources of fluctuations are given by 

C p(sf) (P,r,i) = i^C p(sf) (P,M), (12) 

a 

C p(sf) (p,r,i) = ^< mp(sf) (p,r,t). (13) 

a 

In the following we assume that all the quantities are 
sharply peaked around the Fermi energy and instead of p use 
the quantities e the energy and n the direction of the Fermi 
momentum. Then, for elastic scattering of electrons the fol- 
lowing relation holds 

nwr*"P(*)(p jP ' jr ) = —5(e-e')w imp< - si \n,ri,r), (14) 

where No is the density of states in the Fermi level. 

For a diffusive conductor we apply the standard diffusive 
approximation to the kinetic equations where the charge 
and spin distribution functions are split into the symmetric and 
asymmetric parts: 

/c(s)(n,£,r,t) = / c ( s )o(e,r,f) + n.f c(s)1 (e, r, t). (15) 

Substituting this form of f c in Eqs. and averaging sub- 
sequently over the Fermi momentum direction first weighted 
with one and then with n, we obtain 



V.f cl = J d< f (n,e,r,t), 



Vp 

T 



-g-V/cO — ——fcl 



(16) 



ndn[C P (n, e, r, t) + £ f (n,e, r, t)}. (17) 
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In the same way from Eqs. Jl 5i and 0, for / s , we obtain 

^V.f s i =--/ s0 + / dng(n,e,T,t), (18) 
a TO J 



VF X7f - fsl 
3 3r s 



potential and current density, 

V.Jc - 0, 
V.5j c - if 
jc = -aVipc, 



+ / ndn[C P (n, e, r, t) + £ (n, e, r, t)}, (19) which also imply 



where different relaxation times of normal impurity and spin- 
flip scatterings are defined as 



1 _ 1 1 

111 



(20) 
(21) 



' imp 



= / dnV mp (n,n',r)(n-n'), (22) 



= / dn'w' sf (n,n',r)(n±n'). (23) 



Here we used the identity u; imp(sf )(n, n',r) = w imp(sf) (|n - 
n'|, r), for the elastic scattering. 

In obtaining Eqs. (1161 1 9i we have disregarded terms 
(d/dt)f c t s ) in the expressions (d/dt)f c t s ) — (d/dt + 
WFn.V + eE.V p )/ c ( s ), since we are only interested in the 
zero frequency noise power. The terms of the electric field 
eE. V p are eliminated by substituting e by e — e<^ c ( s ) (r, t) 
in the arguments of / c ( s ), respectively, where the charge and 
spin potentials are expressed as cp c — X^qVq/^ and ip s — 
J2 a aip a /2, with Lp a (r, t) = J dsf a0 {s, r, t), being the spin- 
dependent electro-chemical potential. We also used the iden- 
tities J <in£™ p = 0, which follows from the conservation of 
number of spin a electrons in each normal impurity scattering 
process 9 . In contrast to this, we note that the integral J dn£f^ 
does not vanish, reflecting the fact that spin is not conserved 
by the spin-flip process. 

Combining Eqs. d!6i . (I17> and Jl 8i . dl9l the equations for 
the symmetric parts of the mean charge and spin distribution 
functions are obtained as 



V 2 /co - 0, 

Y72 f _ /s0 
V /so - -T2", 

C sf 



(24) 
(25) 



where 4f — \/D s T s f is the spin-flip length. Note, that in 
general charge and spin diffusion constants given by -D c ( s ) = 
VpT c u\/3 are different. 

Using Eqs. dl5l the charge and spin current densities can be 
expressed as j c ( s ) = (eNgVp/3) J £tef c ( s )i- The corresponding 
fluctuating potentials are given by </? c ( s )(r) + 8ip c ^(r,t) = 
(1/e) f cfe/ c ( s )p. Using these identities and integrating Eqs. 
( I16I17> over e we obtain diffusion equations for the charge 



(26) 
(27) 
(28) 
(29) 



(30) 



In the same way Eqs. dl8H9l > give us diffusion equations of 
spin potential and current density: 



V.Js = (f s 

T s f 

e 2 UF 

T s f 



<Vs + if, 



j s = -er s V<p s , 



V 2 ^ 



9s_ 

t 2 ' 



(31) 



(32) 

(33) 
(34) 

(35) 



where cr c ( s ) = e 2 N D c ^ are the charge and spin conductivi- 
ties. Here 



j c c (s) = ev F N r c{s) J (Q + q s) )ndnde (36) 

are the Langevin sources of fluctuations of the charge and spin 
current densities, and 



if (s) (r,t) = eN J dedn£ {s) 



(37) 



are additional terms in the expressions for the divergence of 
charge and spin currents fluctuations Eqs. d27i and d32l >. due 
to the non-conserved nature of spin-flip process. 

Now we calculate possible correlations between the cur- 
rents , and i% s y From Eqs. d4!6l > for the correlations of 
different fluctuating sources we obtain 



< 



C p(sf) (n,e,r,t)C P(sf) (n',e',r',i') 



>= 



1 



x 6{r - r')6(t - t')6(e ~ e')G™$ sf) (n, n', r, e), (38) 
<C P (n,e,r,i)^(n',e',r',i') >= 0, (39) 



where 



G' Q T = <W J dn"[S(n - n') - 5(n' - n")] 
[J QQ (n, n", e) + J QQ (n", n, e)], (40) 
Gf Q , = J dn"[S aa ,S(n - n') - 6 a _ a ,5(n' - n")] 
[J_ aQ (n,n") + J _ a (n",n)]. (41) 
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v 1= v 



V 2 =0 



FIG. 1 : A schematic picture of the three terminal spin valve. A dif- 
fusive wire of length L is connected to three ferromagnetic terminals 
via the tunnel junctions. 



junction area. Here f a0 = / c o + a/ s o andj Q = j c + cy s are the 
spin a symmetric part of the distribution function and current 
density respectively. From the solutions of the diffusion equa- 
tions the mean charge and spin currents and the correlations 
of the corresponding fluctuations can be obtained. 

Eqs. fl24l35b and J42I47I > are valid for an arbitrary Ti mp /T s f 
in the diffusive limit. In the following we will consider the 
more realistic case of T; mp <C r s f, where the effect of spin- 
flip scattering on the conductivity of the diffusive metal is ne- 
glected. In this case, r s = r c = T; mp and hence er s = er c = a. 
For simplicity, we also assume that the spin-flip scattering is 
isotropic, i. e., w sf does not depend on the directions n, n', 
which implies = r s j" = 2r s f. In the next section we use the 
above developed formalism to calculate spin-polarized current 
correlations in a diffusive three-terminal system. 

III. THREE-TERMINAL SPIN VALVE 



From Eqs. dl 2i . Jl 31 and ( 13814 11 we calculate the corre- 
lations between the fluctuating sources , which can be 
used together with Eqs. j36t and d37l to obtain the results, 

< ie(s)*( r > *)ic(s) m ( r ' ! >= Si m S(r - r')S(t - f>„ w 



x 



Tc(s) 



E[ n -»( r 



' imp 



n„- Q (r)], 



(42) 



< &(r, t)f sm {v', t') >= S lm 6(r - r')6(t - t')a c 
x^an aa (r), (43) 

a 

< i f(r,iK f (s) (r',i')> 

=<Jc(s) m (r , ,O4)( r '' t ')>= . C 44 ) 
< if (r, t)if(r', f ) >= 6 lm 6(r - r')5(t - t')* c (45) 

-^^Tn^r), (46) 



where 



IW(r) = J dsf a0 (s,r)[l - f a >a(e,r)]. 



(47) 



The diffusion equations (I24I35> and Eqs. ( I42I47I > are a com- 
plete set of equations, which in principal can be solved for a 
multi-terminal mesoscopic diffusive conductor connected to 
an arbitrary number N of metallic and/or ferromagnetic ter- 
minals held at constant potentials. The distribution function 
f n (e — eV n ) of the terminal n(— 1, N) biased at the volt- 
age V n determines the boundary conditions of the diffusive 
equations. In the case where the terminal n is connected 
by a tunnel junction to the diffusive conductor at the point 
r„, the ferromagnetic character of the terminal can be mod- 
eled by a spin-dependent conductance g na . The fluctuating 
spin a current through the junction is given by I na (t) = 
g na J de [f a o {e, F n , t) ~ f n [e - eV n )] . As the boundary condi- 
tion this current should match to the value calculated from the 
diffusive equations I na = J A dS ■ j Q (r n , t), where A n is the 



We consider the three-terminal spin valve system as shown 
in Fig. ^ The system consists of a diffusive normal wire (N) 
of length L connected by tunnel junctions to three ferromag- 
netic terminals F^ (i = 1, 2, 3). The terminal Fi is held at the 
voltage V and the voltage at the terminals F2.3 is zero. The 
tunnel junction i connecting Fj to N has a spin dependent tun- 
neling conductances g; Q , which equivalently can be character- 
ized by a total conductance gi = ^ Q gi a , and the polarization 
Pi = a 9ia/gi- Inside the wire we account for both, nor- 
mal impurity and spin flip scattering. The length L is much 
larger than £[ mp providing a diffusive motion of electrons. The 
spin-flip length £ S f is assumed to be much larger than £j mp , but 
arbitrary as compared to L. We study the influence of the spin- 
flip scattering on shot noise of the current through the wire and 
cross-correlations between currents through the terminals F2 
and F3 . 



A. Charge and spin currents fluctuations 



To start we write the solutions of Eqs. J30l > and d35l in 
terms of the charge (spin) potentials <£ c (s)(0) an d <Pc($)(L), at 
the connecting points x = and x = L inside the wire: 

<?c( s )(x) = 0c(s)o(»)^c(s)(O) +<l>c(s)L(%)<Pc(s)( L ), (48) 

where the charge and spin potential functions are defined as 



</>cu(z) = 1 



L 



0sl(z) = 



sinh [A(l - x/L)] 

sinh A 
sinh (Xx/L) 
sinh A 



(49) 
(50) 

(51) 
(52) 



Here the parameter A = L/£ s f is the dimensionless measure 
of the spin-flip scattering inside the N-wire. An expression 
for the fluctuations of the current through the wire A7 c i is 
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obtained if we take the inner product of <5j c in Eq. ( I29l l with 
V(/> c o and integrate over the volume of the wire: 

Aicj =-a J ds- V0 cO <5^ c + 6I£, (53) 

5I C C = J dfi(»f +j c c -V)4o, (54) 

where we also used Eqs. \21\ and J30l >. In a similar way by 
volume integration of the products V0 sO (l) ■ <5js in Eq. ( I34t . 
and using Eqs. d32i and J35i we obtain A/ s i(0, Z), which 
yields the fluctuations of spin currents at x = and x = L: 



the point a; = 0, they read 



AI s i(0, L) = —a J ds ■ V0 8 (o,i)(fy>s 
+« S C (0,L), 

<5/ s c (0,L) = / dO(if +j s c -V)0 sO ,L. 



(55) 
(56) 



Note, that as a result of the spin-flip scattering the spin current 
and, hence, its fluctuations are not conserved through the wire. 
Using Eqs. J55l >. J56b . 15 H and d52l > a relation between the 
fluctuation of the spin currents at the two different points is 
obtained, 



A/,i(0)=AJ,i(i) 

+5i:(o)+5i:(L), 



5J vA 2 t[^ a (0) + ^ B (i)] 



(57) 



where g^ — aA/L (A being the area of the wire) is the 
conductance of the wire and t(A) = tanhA/A. In the limit 
A — > 0, the conservation of the spin current is retained and 
A/ S i (0) = A/ sl (L), as is seen in Eq. J57I . 

In the Boltzmann-Langevin formalism, the fluctuation of 
spin a currents through junction i are written in terms of the 
intrinsic current fluctuations 6Ii a due to the random scattering 
of electrons from the tunnel barriers and the potential fluctua- 
tions 6ip a (0, L) at the junction points: 



A/i a (0, L) = SI ia - gi5ip a {0 7 L). 



(58) 



Using this relation the fluctuations of charge and spin cur- 
rents through the terminals can be expressed in terms of the 
fluctuating spin and charge potentials at the connection points 
and the corresponding intrinsic currents fluctuations. Denot- 
ing SI c r s )i as the intrinsic fluctuations of the charge (spin) cur- 
rent through the tunnel junction i, we obtain 

A7 c i = 5I c i - gidip c (L) - gipiScp s (L), (59) 

Aic2,3 = £J C 2,3 " 32,3<fy?c(0) - 32,3P2,3<Vs(0), (60) 

AJ s i(L) = 5I s i - gipidip c (L) - giSip s (L), (61) 

A/ s2 ,3 = W s 2,3 - 32,3P2,3<5<A:(0) - 32,3<fy>s(0). (62) 

Now we have to apply the currents conservation rules at the 
junction points. For spin-conserving tunnel barriers charge 
and spin currents fluctuations are conserved. At the pint x = 
L the rules apply as equality of the expressions for A/ C i given 
in Eqs. <|53} and (55), and AI sl (L) in Eqs. «|55) and E}. At 



^AI C i = 0, 

i=l 
3 

^A/ Si (0) = 0, 



(63) 



(64) 



which in combination with Eqs. J53l57i and J59l62t lead to 

SI ci = 323^(0) + gi5ip c (L) 



+323P23<fy>s(0) + giPiSip s (L) 

SI»i = 323P23^c(0) + gipiSip c (L) 

i 

+ (323 + 9Nt)6(p s (0) + (gi + g N t)6(p s (L) 
-« S C (0)-« S C (L) 

SIsi = giPiSip c (L) —6<p s (Q) 

s 

8 Id = -gNStp c {0) 

+ (31 + 9n)S(Pc{L) + giPi5ip s (L) - 61% 



(65) 



(66) 



(67) 



(68) 



The solution of this system of equations gives us the fluctua- 
tions of the potentials in the connecting nodes which can be 
replaced into Eqs. ( 1591621 to obtain the fluctuations of the 
charge and spin currents through different terminals in terms 
of SI C i, 5I S i, <5/ c c , and W s c (0, L). In particular, the fluctuations 
of the charge currents have the form 



j'=i 

+c i 6IZ + c i0 6I°(0) + c iL SI°(L), 



(69) 



where c,j, cf.-, Cj, Cjo, and Cj£ are functions of gi, pi, and A. 



B. Mean currents and correlations of currents fluctuations 

The currents correlations (AJ C iAJ C j), are expressed in 
terms of the correlations of different fluctuating currents ap- 
pearing in Eqs. ( I69> . To calculate the correlations of the cur- 
rents 51^ and <5/ s c (0, L) we have to determine the mean distri- 
bution function f a Q = f c o + af s o- This is achieved by solving 
Eqs. (I24i . J25I and imposing the boundary conditions that f a Q 
in the terminal Fj held in equilibrium at the voltage Vi is given 
by the Fermi-Dirac distribution function /$ = Ffd(s — eVi). 
From the solutions of Eqs. (I24> . (125 \ we obtain 

U = fi + (f2-h)[a + b^ 

+a(csinh — + dcosh ^y-)], (70) 

where a, b, c, d are coefficients which have to be determined 
by the boundary conditions. Integrating of d70t over the en- 
ergy e the mean electro-chemical potential of spin a electrons 
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is obtained: 



ip a (x) = [a + b— 

+a(csinh ^ + dcosh -y-)]V, 



(71) 



which also could be obtained from the solutions of Eqs. (130b . 

In the presence of the tunnel junctions the boundary condi- 
tions are imposed by applying the mean currents conservation 
rules at the connection points. Using Eqs. ( I28> . ( I33> . ( I71l > we 
obtain the mean charge and spin currents through the N wire: 



h\ = bg N V, 

Isi{x) =5 , ArA(csinh 



Ax 



dcosh -jr-)V. 
L 



(72) 
(73) 



In terms of the charge and spin potentials at the connection 
points <^ c (0) = aV, (p Q {L) = (a + bL)V, <^ s (0) = dV, 
<p$(L) = (Xsc + cosh Ad) V, we have the following relations 
for the mean currents 

la = -gi<p c (L) - gipi<p s (L), (74) 

Ic2,3 = -ff2,3^c(0) - 52,3P2,3<^s(0), (75) 

hi = -9xPi<Pc{L) - g\(p t {L), (76) 

I S 2,3 = -ff2,3P2,3<^c(0) - 52,3^ S (0), . (77) 

Using Eqs. 1721771 and the currents conservations relations at 
the point x = 0, Yli=i hi = 0, and J>2i=i ^si(O) = 0, we find 

a ,9n . , 9Ns .91 i , , 
77 = [ h (H )?23 H Qi\ cosh A 

L- 523 91 523 



,9N 



9N 



_(l + £il)A 2 
323 91 
9N 

-P1P23- 



(— + qi)q23]s 

9n 



923 
C ff23 

x cosh A - 



(78) 



(1 



9 N \ 1 9l /, 

]923 H (1 

51 323 



323' 



9N_ 

923 



(1 



9n_ 

91 



9n_ 

323' 



)A 2 s 



\qi + 923(1 
-(923 ^ 



b 

C 



9n 
3i 

323 



<7i)]s - 2pip 2 3, 
(71) cosh A 



ran ,2 , ai 1 

-[ — A H <7i<?23R 

323 3w 

c = _9L ( 9_n X 2 s + 

C 323 31 

^ = ^LA^coshA 

L- 323 31 



gi coshA)p 2 3 - 923P1, 
+ gis>23 Api, 

323 

32 + 33,91 = 1 



(79) 

(80) 
(81) 
(82) 



Pi 



where p 2 3 = (32P2 + 33P3V323, 323 
923 = 1 ~P23^ S W = sinhA/A. 

Replacing f a Q given by Eq. J70I in Eqs. J47I we can calcu- 
late the correlations in Eqs. ( I42I46I >. which can be used to cal- 
culate all the possible correlations between 5I£ and <5/ s c (0, L), 



given by Eqs. d54i and (1561 . Calculations lead to the follow- 
ing results 



S c = <« c c « c c ) - 2g N [a(l - a - b) + |(1 - |&) 
H s(caXs H — cosh A)J, 



S.(0) = («J(0)J/f(0)) = 2 9 „[ ? (1 - a) 

c 2 -d 2 1, „ c 2 " 
H — -5 (cdXs H 

5 s (L) = (« s c (L)^(L)>=2. 9Ar [^(l-a) 

c 2 +d 2 

-(cdXs H — cosh A)], 

s 2 



(83) 



■ cosh A)], 



(84) 



2s 2 

S s (0L) = (5i:(0)Si:(L)) = -2g N [- t {\ a) 



(85) 



b . . ,1/t-l 



-d 2 



- —)b 2 + 

2s ! 2st 



1 c 2 + d 2 
- — (cdXs H — cosh A)] , 



(86) 



S(0) = (S^SI^O)} = g N [(l - 2a - 2b)d + (1 - 2a 

cosh A 1 cosh A 

-6)(Ac+ d) + -fec+(d c)6, (87) 

s A As 

S(L) = = -3jv[(1 - 2a - &)(A SC 

l + scoshA A .coshA 1W , coshA ..,.„„. 

+ d) - ( l)(sd+ — : — c)6 .(88) 

s s X 

In writing these equations, we have for simplicity dropped the 
time dependence of the correlators and implicitly assumed, 
that the correlators are symmetrized (which leads to a factor 
of 2). Since we will be solely interested in the zero-frequency 
noise correlations, we also dropped the time integration. 

To obtain the noise in the terminals we have to specify the 
correlators of the intrinsic fluctuations at the tunnel junctions. 
Assuming the tunnel junctions to be spin conserving, we ob- 
tain for the correlations of the intrinsic fluctuation of charge 
and spin currents 



(SI ci SI c j) = (5I si 5I s j) = 
(SI ci SI sj ) = 5ij2e{\I ci+ \ 



(89) 
(90) 



where 7 CJCt = J C j+a/ S j is the mean current of spin a electrons. 

Using the re sults (I83l90l and Eqs. J59I62I with the solu- 
tions of Eqs. ( I65I68> the correlation of the currents of the 
form Sij = (A7 C jA/ CJ -) is obtained. In terms of the coeffi- 
cients introduced in d69i it has the form 



= 2e^2[(c ik c jk + c\ k (? jk )I k + (c ik c s jk 

k=l 

+ c \k c jk){\h+k\ - %-k\)] + (kCjS c + CioCjoSsiO) 
+CiLCjLS s (L) + (CiCj + Ci O Cj)S(0) + {CiC 3 L 

+c iL Cj)S(L) + {cioc-jL + c lL c j0 )S s (0L). (91) 
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FIG. 2: Fano factor F versus the spin flip scattering intensity A = L/l s f for a given magnitude of the polarization in F-terminals but different 
configurations and for different values of the tunnel conductances g = Qi/gN- For the configurations with p2 = pz the corresponding 
cross-correlation is obtained via S23 /] Jci | = (F — l)/4. 



In this way we obtain the Fano factor F = 5n/2e|/ c i| and the 
cross correlations measured between the currents through 
F2, F3. In the general case for arbitrary gi, pi and A the ex- 
pressions of F and S23 are too lengthy to be given here and 
in the next section we will present analytical expressions of F 
and S*23 in some important limits only. 



IV. RESULTS AND DISCUSSION 



A. Shot noise 

Let us start with analyzing of the shot noise. Fig. |2]illus- 
trates the typical behaviour of the shot noise with respect to 
the spin flip scattering intensity and configuration of the po- 
larizations. Here F versus A is plotted for a given magnitude 
of the polarization p in the terminals for the different mag- 
netization configurations. Different columns I-III belong to 
different values of the tunnel contact conductances g. Clearly, 
for a finite A the Fano factor F is changes drastically with the 
relative orientation of the polarizations. For strong spin flip 
scattering, A ^> 1 the Fano factor reduces to 



For simplicity in the following we will consider the junc- 
tions to have the same tunneling conductances gi/gN = 9, 
(i = 1,2,3). If the amplitude of the polarizations {pi} are 
also the same we distinguish the following different configu- 
rations. The two terminals F2 and F3 have either parallel or 
anti-parallel polarizations. We take the signs of P2 and p% to 
be positive in the parallel configuration. In this case there are 
two different configurations depending on whether p\ is posi- 
tive (parallel X.0P2, P3) or negative (anti-parallel top2, P3)- On 
the other hand for the antiparallel configuration of F2 and F3 
the sign of p\ is not essential and the two cases of p\ and —p\ 
are equivalent. Thus there are three independent configura- 
tions of the polarizations corresponding to p-y = P2 = P3 — p, 
—pi = P2 — P3 = p and pi = —p2 = P3 = p. We denote 
these configurations by " + " , " — " , and " 0" , respectively. 

We analyze the dependence of F and F23 = S^/l-fcil on 
the spin flip scattering intensity A for different configurations 
and amplitude of the polarizations. We show how the spin 
flip scattering affects both F and S23 and produces a strong 
dependence on the magnetic configuration of the F terminals. 



F 



N 



1 5 + 6g + Ag 2 + (8/9) g 3 
3 3 + 6 5 + 4.g 2 + (8/9) 5 3 



(92) 



independent of the polarization of the terminals. Eq. J92t 
is the result for a all-normal metal three terminal system 
(pi = P2 = P3 — 0) 9 which reduces to 5/9 and 1 /3 in the lim- 
its of small and large g, respectively. This is expected since 
the strong spin flip intensity destroys the polarization of the 
injected electrons. 

At finite A the curves belonging to different configurations 
differ from each other and the largest difference occurs as A 
approaches zero. In this limit, using Eqs. i9l\ and (1741 we 
obtain the following results for the different configurations 



F, 



T~ 

g X, 



9 



2<? 2 
—9 



16G% 3 , 21. 5 



;)9- 



51 



-q[q{22 + 7q)g + 36 + y q]g 3 



01 94^ 4 

-q[(6q - 32) 5 2 + ll.g + 6] + — (-g 2 + 3g + 1)}, (93) 



8 



F - = ^T[ 10 9 3 ( 2 ~ «).9 3 + V(7 ~ 2q)g 2 + 6g(V 

1 3 RC 1 ^ 2 

-6g + 16) 5 + 3(-g 2 -8q + 16)] + -j^bgV 
4 qx_g z o 
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(94) 
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FIG. 3: Fano factor F versus the magnitude of spin polarization p 
of the terminals in the limit of small spin flip-scattering intensity 
A = L/lsf — > 0. The results are shown for different configurations 
of the polarizations and different values of the tunnel conductances 

g- 



Here we defined x± = 2qg + 3, xq = 2qg + q + 2, and 
q = 1 — p 2 . The total conductances of the system normalized 
by gjv for the three configurations in the limit A — > are given 
by 



G+ 
G_ 
Go 



2gg 2 + 6g + 9/2 ' 
2qgx- 



\2qg + q 

gx 



2qg 2 + 2qg + 4g + q/2 + 4 



In Fig. [3]we show the polarization dependence of the Fano 
factor in the limit of small spin-flip scattering intensity, A — > 
for the different magnetic configurations of the terminals. At 
p = the Fano factor for different conductances g takes the 
normal state value Fn, see Eq. (I92i . independent of the po- 
larizations configuration. For finite polarizations F of differ- 
ent configurations differ from each other and the normal state 
value. As p increases, the Fano factor of the configuration — 
deviates substantially from those of the other two configura- 
tions and reaches the full Poissonian value 1 as p approaches 
1. This is independent of the tunnel conductances g. Thus for 
perfectly polarized terminals the Fano factor takes the value 1 
in the limit of small spin-flip scattering rate. 

To understand this effect we note that in the limit p — > 1, the 
system constitutes an ideal three terminal spin valve due to the 
antiparallel configuration of the polarizations at its two ends. 
In the absence of the spin flip scattering there in no current 
through the N-wire since for the up-spin electrons provided 
by the terminals F 2 and F 3 , there is no empty state in the ter- 
minal Fi in the energy range eV. For very small but finite 



A only those of electrons which undergo spin-flip scattering 
once can carry a small amount of current. These spin-flipped 
electrons are almost uncorrected and pass through the normal 
wire independently giving rise to full Poissonian shot noise. 
Similar effects have been found before for twc 14 and four 18 
terminal spin- valve systems. 

For arbitrary g and p, F has a complicated dependence on 
A and the corresponding expressions are to lengthy to be pre- 
sented here. Simpler expressions are obtained for perfectly 
polarized junctions. Setting p = 1 in Eqs. ( 19 11 1 and J74i yield 
for the Fano factors of the different configurations 



8G^ 



.9 



5G 2 



2 9 



2a 2 

1 + 2(. 9 2 + + 3 5 + 1) 



tanh(A/2) 



cosh A — 3s 
A 2 ^ 2 



3 . tanh 

r a 2 



A 2 

(A/2), 



A 



(96) 



F_ = 



5G 2 



+ - 



9 

8G 3 1 



— Q 

2 y 



l + 2( 3 



2 +2.9 2 
A 2 



3.9+1 



coth (A/2) 



A 



cosh A + 3s 
A^s 2 



3 coth 2 (A/2) 
r A 2 J 



(97) 



9 



X. 



9 
4 

8 Go 



3G ,2 3 



9 



1) 



9 

32 (17-A 2 )+ ff (l-^A 2 ) 9 



A 2 )]- 



(98) 



1 ' " ' ' JX 1 16 

Here we defined x — 2g cosh A + A 2 s(A) and the dimension 
less total conductances are now given by 
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G + = 



G_ 



Go = 



2. 9 + 4 ff tanh(A/2)/A + 3 ' 
1 

2. 9 + 4 3 coth(A/2)/A + 3 ' 



,ga: + 4 5 coshA + (2 ff 2 + 3A 2 /2)s(A) ' 

The strong dependence of F on A, and the magnetic con- 
figuration is also shown in Fig. |4j where F versus A is plot- 
ted for different polarization p\ of the terminal Fi, and fixed 
Pi = P3 — 1- As in Fig. ^different columns I-III present re- 
sults for different values of g. In each column p\ varies from 
top to bottom in the interval —1 to 1. In the limit of large A, 
the Fano factor tends to the the normal state value ( I92i deter- 
mined by the conductance g only. The deviations from this 
normal state value at finite A depend on p% and g. 

For small values of g (column III), with decreasing A from 
large values F first decreases below the normal state value Fn 
and then starts to increase again. Thus, there is a minimum of 
the Fano factor F occurring at a value of A which continu- 
ously decreases from A ~ 1 to A = as p\ increases from — 1 
to 1. Decreasing A further, F increases to a maximum value at 
A = 0. The maximum of F at A = continuously decreases 
with increasing p% and becomes a minimum when p% = 1. 
For antiparallel fully polarized terminals, i. e. p\ — —1, de- 
creasing A leads to the strongest variation of F, see Eq. (197 \ 
and F reaches the Poissonian value 1 as A approaches zero. 
Comparing plots for p\ = —1 of columns (I-III), we see that 
this effect is independent of the contact conductances g, which 
is in agreement with the discussion in connection with Fig. [3] 
Increasing p\ from — 1 decreases the effect of spin flip process 
in the current and the noise, and, hence, leads to a reduction 
of the shot noise at A <C 1. The maximum value of F thus 
drops below the Poissonian value. 

For large values of g (column I) the maximum of F is 
shifted from zero to a finite A ~ 1 as p\ increases from — 1 to 
1, while for small g the maximum F always occurs at A = as 
described above. For g ~ 1 (column II) the situation is in be- 
tween the large and small g behavior, increasing A leads first to 
a maximum of the Fano factor followed by a minimum. Com- 
paring the first and the second row of Fig.0] we observe the 
effect of spin-flip scattering is more pronounced for negative 
pi than for positive p\. This can be understood, because for 
parallel magnetization directions the non-equilibrium spin ac- 
cumulation in the ferromagnetic wire is smaller. The spin-flip 
scattering decreases the spin-accumulation and, hence, has the 
largest effect for anti-parallel magnetizations. 



B. Cross-correlations 

Let us now discuss the effect of spin flip scattering on the 
cross correlations measured between the currents through the 
terminals F2 and F3. We distinguish two different cases of 
parallel (p2 = P3) and antiparallel (p2 = — P3) magnetizations 
of F2 and F3. For the parallel case and when 52 = 33 the two 
terminals F2, F3 are completely equivalent and hence AI2 = 
A/3. This can be used with Eq. d63i to obtain that in this 
case the Fano factor and the cross-correlation factor F23 = 
523/|2e/ c i| are related as (see also Refi2£) 



F23 — 



F- 1 



(99) 



Thus F23 has the same qualitative dependence on A as F. 
Since F < 1 the cross-correlations are always negative as 
expected 31 . The Fano factor for the perfectly polarized par- 
allel case is presented in Fig.|4]and the cross-correlations for 
this case can be deduced from these plots using Eq. i99l . We 
will now analyze the cross-correlations F23 for P2 = P3 = 1 
and different values of p\. From Eqs. ( 19 U and J74i one can 
see that in the limit of large spin-flip scattering A ^> 1, the 
cross-correlations reduces to its all-normal system value 



F23 — — — 



l (4/3)g 3 + 6g 2 + 9g + 3 
9(8/9) 5 3 + 4 3 2 + 6 ff + 3' 



(100) 



which is independent of the polarizations. Alternatively, this 
result could have been obtained using Eqs. 03 and J9§. On 
lowering A the amplitude of the cross-correlations I-F23I de- 
creases (large g, column I) or increases (small g, column III) 
with respect to the normal value. For p% = — 1, F23 vanishes 
in the limit A — > 0, irrespective of the the value of the contacts 
conductance g. This case corresponds to the vanishing of the 
total mean current I\. From this observation we conclude that 
in the expression of F23, the cross-correlation S23 vanishes 
faster than the mean current I\, as A — > 0. Thus, for p\ = — 1 
I-F23I has a minimum at A = for both cases of small and 
large g. Increasing p\ from —1, the minimum is shifted to a 
finite A ~ 1 for large g, it stays always at A = for small 
g. For small g, |i<23| has also a maximum at A ~ 1 which 
corresponds to the minimum of the Fano factor. 

For g ~ 1 (column II) the behavior of F23 is between that 
for small and large g. Comparing of the plots in top (pi < 0) 
and down (pi > 0) rows in Fig. |4]shows that the most strong 
variation of the cross-correlation happens when the magneti- 
zation of Fi is anti-parallel to those of F2,3. 

For the antiparallel case P2 = —P3 the effect of spin-flip 
scattering is more interesting, since it contains the correlations 
between currents of opposite spin directions produced by the 
spin-flip scattering. In Fig. [5] we plot F23 versus A for different 
values of the magnitude = \ps\ — P and g. We take 
pi = 1 which corresponds to a maximum spin accumulation 
in the N-wire due to the terminal Fi . For p = 1 the cross- 
correlations are solely due to the spin flip-scattering. In this 
case ^23 vanishes in the limit of A — > 0. 

At finite A the spin-flip scattering induces correlations be- 
tween the electrons with opposite spins and hence F23 be- 
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FIG. 4: Fano factor F versus the spin flip scattering intensity A = L/£ s f for different polarizations pi in the terminal Fi, and perfect 
polarizations in the other two terminals F2, F3: P2 = pz = 1. The columns I-III correspond to different values of the tunnel contact 
conductances g = gi /gjv. Here the cross-correlations are related to the Fano factor via SizlXlti | = (F — l)/4. For explanation of the various 
plots, see text. 



comes finite. With increasing A the magnitude of F23 in- 
creases and approaches the all-normal system value of pi = 
P2 = Pz = when A > 1. For p < 1 the value of |F 2 3| for 
vanishing A depends on the values of p and g as 



F 23 = --(l-p 2 



(4/3)<? 3 + 4<? 2 + 4.g+l 



(101) 



(.9 + I) 3 

Decreasing p from 1 to 0, |F 2 3 1 increases from zero to a max- 
imum value. The maximum absolute value is equal to the 
normal value for large g (column I), while it is larger than the 
normal value for small g (column III). 



V. CONCLUSIONS 

In conclusion, we have investigated the influence of spin 
polarization and spin-flip scattering on current fluctuations 
in a three-terminal spin- valve system. Based on a spin- 
dependent Boltzmann-Langevin formalism, which accounts 



for spin-flip scattering in addition to the usual scattering at im- 
purities and tunnel junctions, we have developed a semiclassi- 
cal theory of current fluctuations in diffusive spin-valves. This 
theory allows the calculations of spin-polarized mean currents 
and correlations of the corresponding fluctuations in multi- 
terminal systems of diffusive wires and tunnel contacts. 

We have applied this formalism to a three-terminal system 
consisting of a diffusive normal wire connected at the ends to 
one and two ferromagnetic terminals, respectively. We have 
found a strong deviation of the current correlations from the 
all-normal system values. The shot noise of the total current 
through the system and the cross-correlations between cur- 
rents of two different terminals depend strongly on the spin- 
flip scattering rate and the spin polarization and change dras- 
tically with reversing the polarizations in one or more of the 
terminals. 

The strongest variation of the shot noise occurs, when the 
polarizations of the two terminals connected to one end of the 
normal wire are antiparallel with respect to the terminal on the 
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FIG. 5: Cross-correlations of the current measured between terminals F2 and F3 for anti-parallel configuration of the polarizations p2 = —pa 
and for different contact conductances g — Qi/gN- The polarization of Fi, pi = 1 and the magnitude of the magnetizations |pa| = | j03 1 vary 
in each column. 



other end. For small spin-flip scattering intensity, the Fano 
factor deviates substantially from the normal value and can 
reach the full Poissonian value for perfectly polarized termi- 
nals even if the tunnel contact resistances are negligible. 

We have further demonstrated the effect of spin- 
polarization and spin-flip scattering on the cross-correlations 
measured between currents of two adjacent terminals in two 
cases where the terminals have parallel and antiparallel polar- 
izations. For antiparallel orientations of the contact polariza- 
tions, the noise allows a direct determination of the spin-flip 
scattering processes in the normal wire. The study of noise 
and cross-correlations therefore allows to extract information 



on the spin-flip scattering strength, which is of importance for 
spintronics applications. 
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